Comparison of two model frameworks for fiber dispersion in the elasticity of soft biological tissues

This study compares two models that are used to describe the elastic properties of fiber-reinforced materials with dispersed fibers, in particular some soft biological tissues such as arterial walls and cartilages. The two model approaches involve different constitutive frameworks, one being based on a generalized structure tensor (GST) and the other on the method of angular integration (AI). By using two representative examples, with the same number of parameters for each model, it is shown that the predictions of the two models are virtually identical for a significant range of large deformations, which contradicts conclusions contained in several papers that are based on faulty analysis. Additionally, each of the models is fitted to sets of uniaxial data from the circumferential and axial directions of the adventitia of a human aorta, both models providing excellent agreement with the data. While the predictions of the two models are comparable and exclusion of compressed fibers can be accommodated by either model, it is well known that the AI model requires more computational time than the GST model when used within a finite element environment, in particular if compressed fibers are excluded

Modeling of fibrous biological tissues with a general invariant that excludes compressed fibers

Dispersed collagen fibers in fibrous soft biological tissues have a significant effect on the overall mechanical behavior of the tissues. Constitutive modeling of the detailed structure obtained by using advanced imaging modalities has been investigated extensively in the last decade. In particular, our group has previously proposed a fiber dispersion model based on a generalized structure tensor. However, the fiber tension–compression switch described in that study is unable to exclude compressed fibers within a dispersion and the model requires modification so as to avoid some unphysical effects. In a recent paper we have proposed a method which avoids such problems, but in this present study we introduce an alternative approach by using a new general invariant that only depends on the fibers under tension so that compressed fibers within a dispersion do not contribute to the strain-energy function. We then provide expressions for the associated Cauchy stress and elasticity tensors in a decoupled form. We have also implemented the proposed model in a finite element analysis program and illustrated the implementation with three representative examples: simple tension and compression, simple shear, and unconfined compression on articular cartilage. We have obtained very good agreement with the analytical solutions that are available for the first two examples. The third example shows the efficacy of the fibrous tissue model in a larger scale simulation. For comparison we also provide results for the three examples with the compressed fibers included, and the results are completely different. If the distribution of collagen fibers is such that it is appropriate to exclude compressed fibers then such a model should be adopted

Classical and all-floating FETI methods for the simulation of arterial tissues

High-resolution and anatomically realistic computer models of biological soft tissues play a significant role in the understanding of the function of cardiovascular components in health and disease. However, the computational effort to handle fine grids to resolve the geometries as well as sophisticated tissue models is very challenging. One possibility to derive a strongly scalable parallel solution algorithm is to consider finite element tearing and interconnecting (FETI) methods. In this study we propose and investigate the application of FETI methods to simulate the elastic behavior of biological soft tissues. As one particular example we choose the artery which is - as most other biological tissues - characterized by anisotropic and nonlinear material properties. We compare two specific approaches of FETI methods, classical and all-floating, and investigate the numerical behavior of different preconditioning techniques. In comparison to classical FETI, the all-floating approach has not only advantages concerning the implementation but in many cases also concerning the convergence of the global iterative solution method. This behavior is illustrated with numerical examples. We present results of linear elastic simulations to show convergence rates, as expected from the theory, and results from the more sophisticated nonlinear case where we apply a well-known anisotropic model to the realistic geometry of an artery. Although the FETI methods have a great applicability on artery simulations we will also discuss some limitations concerning the dependence on material parameters.Comment: 29 page

Modeling and experimental investigations of the stress-softening behavior of soft collagenous tissues

This paper deals with the formulation of a micro-mechanically based dam-age model for soft collagenous tissues. The model is motivated by (i) a sliding filament model proposed in the literature [1] and (ii) by experimental observations from electron microscopy (EM) images of human abdominal aorta specimens, see [2]. Specifically, we derive a continuum damage model that takes into account statistically distributed pro- teoglycan (PG) bridges. The damage model is embedded into the constitutive framework proposed by Balzani et al. [3] and adjusted to cyclic uniaxial tension tests of a hu- man carotid artery. Furthermore, the resulting damage distribution of the model after a circumferential overstretch of a simplified arterial section is analyzed in a finite element calculation

An arterial constitutive model accounting for collagen content and cross-linking

It is apparent from the literature that the density of cross-links in collagenous tissue has a stiffening effect on the mechanical response of the tissue. This paper represents an initial attempt to characterize this effect on the elastic response, specifically in respect of arterial tissue. Two approaches are presented. First, a simple phenomenological continuum model with a cross-link-dependent stiffness is considered, and the influence of the cross-link density on the response in uniaxial tension is illustrated. In the second approach, a 3D model is developed that accounts for the relative orientation and stiffness of (two families of) collagen fibers and cross-links and their coupling using an invariant-based strain-energy function. This is also illustrated for uniaxial tension, and the influence of different cross-link arrangements and material parameters is detailed. Specialization of the model for plane strain is then used to show the effect of the cross-link orientation (relative to the fibers) and cross-link density on the shear stress versus the amount of shear deformation response. The elasticity tensor for the general (3D) case is provided with a view to subsequent finite element implementation

Biomechanical relevance of the microstructure in artery walls with a focus on passive and active components

The microstructure of arteries, consisting, in particular, of collagen, elastin and vascular smooth muscle cells, plays a very significant role in their biomechanical response during a cardiac cycle. In this paper we highlight the microstructure and the contributions of each of its components to the overall mechanical behavior. We also describe the changes of the microstructure which occur as a result of abdominal aortic aneurysms and disease such as atherosclerosis. We also focus on how the passive and active constituents are incorporated into a mathematical model without going into detail of the mathematical formulation. We conclude by mentioning open problems towards a better characterization of the biomechanical aspects of arteries that will be beneficial for a better understanding of cardiovascular pathophysiology

Modeling of damage in soft biological tissues and application to arterial walls

A new material model is proposed for the description of stress-softening observed in cyclic tension tests performed on soft biological tissues. The modeling framework is based on the concept of internal variables introducing a scalar-valued variable for the representation of fiber damage. Remanent strains in fiber direction can be represented as a result of microscopic damage of the fiber crosslinks. Particular internal variables are defined able to capture the nature of soft biological tissues that no damage occurs in the physiological loading domain. A specific model is adjusted to experimental data taking into account the supra-physiological loading regime. For the description of the physiological domain polyconvex functions are used which also take into account fiber dispersion in a phenomenological approach. The applicability of the model in numerical simulations is shown by a representative example where the damage distribution in an arterial cross-section is analyzed

Modeling of damage in soft biological tissues and application to arterial walls

A new material model is proposed for the description of stress-softening observed in cyclic tension tests performed on soft biological tissues. The modeling framework is based on the concept of internal variables introducing a scalar-valued variable for the representation of fiber damage. Remanent strains in fiber direction can be represented as a result of microscopic damage of the fiber crosslinks. Particular internal variables are defined able to capture the nature of soft biological tissues that no damage occurs in the physiological loading domain. A specific model is adjusted to experimental data taking into account the supra-physiological loading regime. For the description of the physiological domain polyconvex functions are used which also take into account fiber dispersion in a phenomenological approach. The applicability of the model in numerical simulations is shown by a representative example where the damage distribution in an arterial cross-section is analyzed

An affine continuum mechanical model for cross-linked F-actin networks with compliant linker proteins

Cross-linked actin networks are important building blocks of the cytoskeleton. In order to gain deeper insight into the interpretation of experimental data on actin networks, adequate models are required. In this paper we introduce an affine constitutive network model for cross-linked F-actin networks based on nonlinear continuum mechanics, and specialize it in order to reproduce the experimental behavior of in vitro reconstituted model networks. The model is based on the elastic properties of single filaments embedded in an isotropic matrix such that the overall properties of the composite are described by a free-energy function. In particular, we are able to obtain the experimentally determined shear and normal stress responses of cross-linked actin networks typically observed in rheometer tests. In the present study an extensive analysis is performed by applying the proposed model network to a simple shear deformation. The single filament model is then extended by incorporating the compliance of cross-linker proteins and further extended by including viscoelasticity. All that is needed for the finite element implementation is the constitutive model for the filaments, the linkers and the matrix, and the associated elasticity tensor in either the Lagrangian or Eulerian formulation. The model facilitates parameter studies of experimental setups such as micropipette aspiration experiments and we present such studies to illustrate the efficacy of this modeling approach

An exponential constitutive model excluding fibers under compression: application to extension-inflation of a residually stressed carotid artery

Detailed information on the three-dimensional dispersion of collagen fibres within layers of healthy and diseased soft biological tissues has been reported recently. Previously we have proposed a constitutive model for soft fibrous solids based on the angular integration approach which allows the exclusion of any compressed collagen fibre within the dispersion. In addition, a computational implementation of that model in a general purpose finite element program has been investigated and verified with the standard fibre-reinforcing model for fibre contributions. In this study, we develop the proposed fibre dispersion model further using an exponential form of the strain-energy function for the fibre contributions. The finite element implementation of this model with a rotationally symmetrical dispersion of fibres is also presented. This includes explicit expressions for the stress and elasticity tensors. The performance and implementation of the new model are demonstrated by means of a uniaxial extension test, a simple shear test, and an extension–inflation simulation of a residually stressed carotid artery segment. In each example we have obtained good agreement between the finite element solution and the analytical or experimental results